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ABSTRACT 


This  report  analyzes  the  effects  of  local  water  move¬ 
ment  on  the  output  of  sensors  that  measure  ship's  speed 
relative  to  water.  The  four  water  disturbances  considered 
are  surface  waves,  internal  waves,  turbulence,  and  vertical 
gradients  of  currents.  The  effects  range  from  quasi- 
periodic  variations  in  the  indicated  speed  caused  by  wave 
motions  to  offset  errors  caused  by  current  changes  with 
depth.  The  quasi -periodic  variations  cause  difficulty  in 
averaging  and  sampling  the  water  speed  for  use  in  other 
ship's  systems.  It  is  concluded  that  local  water  movement 
does  limit  the  useful  accuracy  of  ship's  water  speed  sensors. 
Evaluating  the  magnitude  of  this  limit  for  any  particular 
ship  and  sensor  and  determining  the  controlling  factor  will 
require  careful  analysis  of  the  specific  case  and  at-sea 
experimentation  to  confirm  the  results. 
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REPORT  HIGHLIGHTS 


OBJECTIVE 

The  objective  of  this  study  is  to  analyze  the  effects  of  local  water  motion  on 
the  useful  accuracy  of  sensors  that  measure  ship's  speed  relative  to  water.  The 
water  motions  considered  arise  from  surface  waves,  internal  waves,  turbulence,  and 
vertical  gradients  of  horizontal  currents.  The  oscillatory  and  offset  errors  in  speed 
signals  affect  fire-control  solutions  and  navigation. 

APPROACH 

The  report:  (a)  describes  the  phenomena  with  mathematical  relations  between 
local  water  motion  (such  as  wave  orbital  velocities)  and  environmental  parameters 
(such  as  wind),  (b)  estimates  the  types  and  magnitudes  of  speed  sensing  errors 
induced  by  the  environments,  and  (c)  presents  ranges  of  errors  a  naval  vessel  could 
expect. 

RESULTS 

The  study  identifies  several  problem  areas  in  speed  sensing: 

•  Surface  waves  cause  oscillating  horizontal  currents  with  peak  magni¬ 
tudes  as  high  as  several  knots.  The  magnitudes  decrease  with  depth,  with  the 
spectral  peaks  tending  toward  lower  frequencies  (or  longer  wavelengths)  at  increasing 
depths.  This  spectral  change  may  cause  difficulties  for  submarines  in  selecting  a 
data  sampling  and  processing  or  averaging  scheme  that  can  reduce  wave  effects  and 
still  respond  to  ship  maneuvers. 

•  The  horizontal  components  of  oscillating  currents  associated  with 
internal  waves  are  proportional  to  the  wave  amplitude  and  phase  speed  and  have  magni¬ 
tude  up  to  tenths  of  knots.  Since  the  wave  amplitudes  do  not  exceed  the  upper  layer 
thickness  (depth  to  thermocline)  and  the  wave  phase  speed  can  be  calculated  or  mea¬ 
sured,  it  should  be  possible  to  estimate  the  upper  limit  to  the  horizontal  currents 
associated  with  internal  waves  by  using  local  measurements. 

•  For  vessels  of  larger  mass  in  relatively  rapid  local  water  movements, 
the  water  motion  causes  fluctuations  in  the  indicated  water  speed  while  the  speed  over 
ground  is  essentially  constant.  For  vessels  of  lesser  mass  in  slowly  varying  local 
water  movements,  the  speed  relative  to  water  is  very  nearly  constant  while  the  speed 
over  ground  varies  as  the  water  motion.  These  results  follow  directly  from  the 
inertial  and  drag  characteristics  of  the  vessel. 
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•  Turbulence  causes  random  but  small  (~0. 01  to  0. 1  knot)  fluctuations 
in  the  speed  readout.  The  spectral  form  appropriate  to  homogeneous,  isotropic 
turbulence  may  be  useful  in  predicting  speed  readout  variations. 

•  Local  ocean  currents  tend  to  decrease  in  magnitude  and  shift  direc¬ 
tions  as  a  function  of  depth.  A  submersible  and  a  surface  support  vessel  operating 
together  will  almost  certainly  experience  differences  between  the  indicated  water 
speeds  of  the  two  vessels  even  while  they  maintain  zero  relative  velocity.  Typical 
gradients  for  the  North  Atlantic  show  about  0.4-knot  change  from  the  surface  to  a 
depth  of  50  meters. 
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NOTATION 


A  -  Empirical  constant  in  equation  (5),  A  =  8.29  x  10“7 
a  -  Wave  amplitude,  equation  (8) 

B  -  Empirical  constant  in  equation  (5),  B  =  4.76  x  10"4 
c  -  Internal  wave  speed  defined  by  equation  (7) 

Cjj  -  Drag  coefficient,  equation  (13) 

D(r)  -  Velocity  structure  function,  equation  (19) 

d  -  Depth  of  water 

E(k)  -  Turbulence  spectra 

f  -  Frequency 

g  -  Acceleration  of  gravity 

h'  -  Mean  depth  of  upper  layer  (see  figure  4) 

h"  -  Mean  depth  of  lower  layer  (see  figure  4) 

k  -  Wave  number  (k  =  2n/\,  where  X  =  wavelength) 

M  -  Mass  of  submarine 
m  -  Number  defined  by  equation  (2)  or  (3). 

N  -  Brunt -Vaisala  frequency,  equation  (10) 

r  -  Separation  distance,  equation  (19) 

S(f)  -  Surface  wave  spectrum,  defined  by  equation  (5). 

T  -  Thrust  of  submarine  propeller. 

T0  -  Mean  thrust 

t  -  Time 

U  -  Horizontal  current  caused  by  internal  waves,  equation  (9) 

Preceding  page  blank 
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tr  -  Subsurface  velocity  variance  caused  by  wind  waves,  equation  (1) 
V  -  Speed  of  submarine  relative  to  ground 
VQ  -  Mean  speed  of  submarine  relative  to  ground. 

Vj  -  Indicated  speed  relative  to  water 
AV  -  Small  perturbation  from  mean  speed 
v  -  instantaneous  current,  equation  (12) 
v0  -  Maximum  amplitude  of  current,  equation  (12) 

W  -  Wind  speed 

z  -  Vertical  coordinate 

e  -  Energy  dissipation  rate 

•  * 

7)  -  Density  discontinuity  surface,  (figure  4) 

9  -  Latitude 

k  -  Turbulence  wave  number 

p  -  Density 

p'  -  Density  of  upper  layer  (see  figure  4) 

p"  -  Density  of  lower  layer  (see  figure  4) 

<r  -  Frequency  of  internal  wave,  equation  (8) 

0  -  Phase  angle  defined  by  equation  (.18) 

u)  -  Radian  frequency  seen  by  observer,  equation  (12) 


L! 


27-446 


x 


i 


INTRODUCTION 


Naval  vehicles  require  a  self-contained  speed-measurement  capability  for 
inputs  to  systems  such  as  the  ship's  inertial  navigator  or  the  fire-control  computers 
as  well  as  for  navigation  when  external  position  fixes  are  not  available.  Most  present 
day  vehicles  have  log  systems  that  measure  the  speed  of  the  ship  relative  to  the  water. 
Thus,  motion  of  the  local  water  in  the  sensing  region  will  cause  a  variance  (error) 
in  the  speed  signal. 

The  purpose  of  this  study  is  to  analyze  the  effects  of  local  water  motion  on  the 
useful  accuracy  of  sensors  that  measure  ship's  speed  relative  to  water.  Four  factors 
are  considered.  Three  are:  surface  waves,  internal  waves,  and  turbulence,  all  of 
which  cause  oscillatory  or  random  errors.  The  fourth  factor  is  vertical  current 
gradients  that  cause  errors  depending  on  the  depth  of  sensing  region.  These  factors 
are  considered  local  in  the  sense  that  they  occur  over  distance  scales  from  tens  of 
feet  to  several  miles  and  time  scales  of  minutes  to  several  hours. 

The  basic  theory  describing  the  water  motions  is  given  in  the  summary  by 
Williams1  with  example  data  taken  from  standard  oceanographic  atlases.  3 

UNCERTAINTIES  CAUSED  BY  SURFACE  WAVES 

The  water  particles  below  a  passing  surface  wave  undergo  an  elliptic  motion 
with  the  size  of  the  ellipse  decreasing  exponentially  with  depth.4  Ideally,  the 
elliptic  path  closes  on  itself  and  no  net  movement  of  the  water  occurs.  In  reality, 
nonlinearties  of  the  wave  motion  exist,  and  there  is  a  small  net  transport  current  in 
the  direction  of  the  wave  movement.  However,  even  with  little  or  no  net  transport, 
the  oscillatory  motion  of  the  water  will  introduce  a  fluctuating  velocity  signal  with 
a  peak  amplitude  that  may  be  as  high  as  several  knots. 4  With  wave  periods  as  long 
as  20  seconds,  these  velocity  fluctuations  may  cause  appreciable  errors  in  computa¬ 
tions  based  on  water  speed.  When  the  vessel  is  moving  in  the  direction  of  the  most 
significant  wave  field,  the  doppler  shift  may  cause  the  apparent  period  to  be  much 
longer  with  the  attendant  difficulty  of  averaging  out  the  wave-induced  water  velocities. 

The  quasi -periodic  oscillations  in  the  indicated  speed  become  very  significant 
when  the  speed  data  are  sampled  at  some  relatively  low  frequency  for  input  into  a 
computer.  If  the  sampling  period  is  harmonic  with  one  of  the  dominant  wave  periods, 
the  computer  averaging  may  be  incorrect.  The  sampling  problem  is  complicated  for 
submarines  because  the  dominant  wave  period  increases  with  depth.  A  sampling 
frequency  which  is  optimum  at  one  depth  is  not  necessarily  optimum  at  another.  A 
compensating  factor,  however,  is  the  exponential  decrease  of  the  wave  currents 
with  depth. 


Superscripts  refer  to  similarly  numbered  entries  in  the  Technical  References 
at  the  end  of  the  text. 
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The  variance  of  the  subsurface  velocities  caused  by  surface  waves  is  calculated 
as  follows. 5  Assume  the  horizontal  velocity  components  have  a  normal  distribution 
with  a  zero  mean.  The  normal  distribution  is  completely  specified  by_its  variance. 
From  hydrodynamic  potential  theory  the  subsurface  velocity  variance  u2  is  related 
to  the  surface  wave  spectrum  S(f)  by 


n 


2irt  cosh  m(d-t-z) 
sinh  md 


2 


S(f)  df 


where  m  is  defined  by 


(2nf)^  -•  mg  tanh  md 


(2) 


f  =  frequency  of  surface  or  wave  spectra,  Hz* 
d  =  depth  of  water  to  bottom 

z  =  distance  from  surface  to  point  u2  is  considered  (z<0) 
g  =  acceleration  of  gravity. 


If  the  water  depth  is  much  greater  than  the  longest  wavelength  considered, 
then  tanh  md  is  very  nearly  unity  and  we  may  write 


m 


47T2f2 

t 

g 


(3) 


Using  the  approximations  for  deep  water  where  d  »  \  ,  the  expression  for  is 
simplified  to 


/ 


47r2f2  exp  —  z]  S(f)  df. 


Z) 


(4) 


The  surface  or  wave  spectra  can  be  described  by  the  Pierson-Moskowitz  spectrum 
for  a  fully  developed  sea.  6»^ 


S(f)  df  =  Ag2f  "5  exp 


[-B(g/fW)4] 


df 


(5) 


♦Abbreviations  used  in  this  text  are  from  the  GPO  Style  Manual,  1973,  Unless 
otherwise  noted. 
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where: 


A  =  8.29  X  10-7 
B  =  4.76  X  10-4 
g  =  acceleration  of  gravity 
W  =  wind  velocity 
f  =  frequency. 

The  variance  of  the  horizontal  component  of  the  subsurface  velocities  caused 
by  surface  waves  is  related  to  wind  velocity  and  depth  by  the  integral  equation: 

u2  =  n 47r2Ag2  exp  [(8?r2  f 2  z/g)  -  (Bg4/f4W4)]  f“3  df.  (6) 

Jo 

Figure  1  shows  sections  of  analog  records  of  the  variations  in  the  speed 
signal  from  an  electromagnetic  speed  sensor  on  a  submarine  at  various  depths 
beneath  a  state  6  sea.  The  curves  shown  in  figure  2  represent  calculations  based 
on  equation  (6).  The  circled  points  were  obtained  from  the  submarine  data. 


UNCERTAINTIES  CAUSED  BY  INTERNAL  WAVES 

The  gravity  waves  observed  on  the  sea  surface  are  a  special  example  of  the 
more  general  case  of  waves  that  may  exist  wherever  there  is  a  vertical  density 
gradient.  Such  a  density  gradient  exists  at  all  times  throughout  an  ocean-water 
column,  but  it  is  greatest  in  the  thermocline.  The  density  gradient  associated  with 
a  thermocline  is  illustrated  by  the  two  examples  taken  from  standard  atlases3  and 
shown  in  figure  3. 


To  simplify  the  mathematics,  consider  a  two-layer  model  where  the  density 
gradient  is  a  step  change  between  the  layers  (see  figure  4).  For  a  shallow  upper 
layer  and  deep  lower  layer  (k  h"  large,  k  h*  small),  the  wave  speed  is  given  by 


°  -  M^)] 


1/2 


(?) 


where 


c  =  wave  phase  velocity 
g  =  gravity 

h'  =  thickness  of  upper  layer 
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p'  =  density  of  upper  layer 


p"  =  density  of  lower  layer. 

Assuming  the  vertical  displacement  of  the  interface  is  a  simple  cosine, 


t)  =  a  cos  (kx-fft), 


where: 


a  =  wave  amplitude 
k  =  wave  number  (k  =  2tt/A) 

<r  =  frequency; 

the  horizontal  flow  in  the  upper  layer  is  given  by 


U  =  -  —  C  cos(kx-at). 


The  frequency  of  the  possible  natural  internal  waves  is  bounded  by  the  Brunt - 
Vaisala  or  stability  frequency  on  the  highest  end  and  the  inertial  frequency  an  the  low- 
frequency  side.  These  frequencies  are 

f  g  dpi1/2  r _ _  .....  , _ 


N  “  [-  P  f  ]  [Brunt-Vaisaia  frequency] 


where 


g  =  gravity 


=  density 


=  mean  density  gradient 


.  _  sin  0 

inertial  12  hours 


[inertial  frequency] 


where  6  =  latitude. 
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The  frequency  and  propagation  velocity  of  the  internal  wave  are  determined  to  a 
large  extent  by  the  density  gradient.  However,  the  wave  amplitude  (hence,  horizontal 
currents)  are  determined  primarily  by  the  generating  source.  The  induced  currents 
are  directly  proportional  to  the  ratio  of  wave  amplitude  to  the  upper  layer  thickness 
as  shown  in  equation  (9).  Since  the  amplitudes  of  the  internal  waves  in  the  ocean  do 
not  exceed  the  upper  layer  thickness  (a<  h'),  the  wave-induced  currents  will  be  less 
than  the  phase  velocity  of  the  wave, 


U  <  c. 

Representative  current  magnitudes  attributed  to  internal  waves  are  of  the  order  of  1 
to  20  cm/s. 8 

The  effect  of  the  oscillating  horizontal  current  of  an  internal  wave  on  submarine 
motion  is  derived  in  appendix  A.  The  drag  forces  of  a  submarine  depend  on  its 
velocity  relative  to  the  water  whereas  the  acceleration  J'orces  depend  on  changes  in 
speed  relative  to  ground.  The  oscillating  currents  of  internal  waves  affect  both  speed 
over  ground  and  speed  relative  to  water. 

Consider  an  oscillating  horizontal  current  component  v,  seen  by  an  observer, 
being  translated  at  steady  speed  Vc.  Thus, 

v  =  v0  sin  cut  (12) 


where 

v  =  instantaneous  current  seen  by  observer 
vQ  =  Maximum  amplitude  of  current 
cu  =  radian  frequency  seen  by  observer 
t  -  time. 

A  submarine  of  mass  M  moving  through  this  flow  at  constant  thrust  T0  will 
undergo  accelerations  because  of  drag  variations  from  changing  flow.  Let  V0  be 
the  submarine  water  speed  corresponding  to  T  ::  T0  with  vQ  =  0.  Assuming  the 
drag  is  proportional  to  the  square  of  the  velocity  through  the  water,  we  have 

M  ~  +  Cd(V-v)2  +  T0  =  0  (13) 


where: 


V  =  speed  of  submarine  relative  to  ground 
Cd  =  drag  coefficient 
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M  =  mass  of  submarine 
T0  =  constant  propeller  thrust. 

By  expressing  the  velocity  of  the  submarine  over  the  ground  as 

V  =  V0  +  AV  (14) 

(where  AV  is  a  small  perturbation  from  the  average  speed  V0)  and  assuming  that 

V0  >>  |v-AV[ 

the  solution  to  equation  (13)  is  found  in  appendix  A  to  be 

V  =  V0  +  v0  cos  V  sin(cot  +  <p)  (15) 


where 

oiM 

tan?  =  -jv^  '  (16) 

The  indicated  speed,  Vj  (or  speed  relative  to  the  water)  is  just  V-v  or 

Vj  =  V0  +  v0  sin  <P  cos(cot  +  <p)  .  (17) 

Using  these  equations  and  the  known  mass  and  drag  characteristics  ol  the  sub¬ 
marine,  the  effects  of  current  fluctuations  on  speed  sensing  can  be  determined.  First, 
the  differences  between  VQ  and  both  V  and  Vj  are  proportional  to  the  peak  current 
magnitude  v0.  Second,  both  V  and  Vj  depend  on  the  angle  <P,  which  is  related  to  the 
submersible  and  current  parameters  through  equation  (16).  As  <P-*o,  Vj-*V0  and 
V  —  VQ  +  vQsin  cot.  The  indicated  speed  approaches  the  constant  valve  V0,  but  the 
speed  over  the  ground  varies  as  the  oscillatory  current.  This  case  represents  a 
small  submersible  with  a  relatively  high  drag  in  a  slowly  changing  current.  The 
other  extreme,  as  <p  —  90°,  results  in  Vj  —  V0  -  v0  sin  cot  and  V  -*  VQ.  This  case 
represents  a  large  submersible  with  low  drag  in  a  rapidly  varying  current. 

The  results  of  an  analog  computer  solution  of  equation  (13)  are  given  in  figure 
5  for  the  assumed  parameters: 

M  =  5,000  tons 

V0  -  5  kn 

T0  =  22, 000  lb 
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0.5  kn* 


u  =  0.0105  s"l  (for  item  (a),  figure  5,  a  10-min  wave  period) 

«  0.00105  s“l  (for  item  (b),  figure  5,  a  100-min  wave  period) 

UNCERTAINTIES  CAUSED  BY  TURBULENCE 

Turbulence  is  an  irregular  fluctuation  of  flow  superimposed  on  the  general  or 
mean  flow  pattern.  Small-scale  turbulence  covering  a  few  tens  of  meters  will 
generate  irregular  fluctuations  in  the  sp^ed  measurement.  While  these  small-scale 
turbulent  eddies  cause  speed  fluctuations  that  tend  to  average  to  zero,  there  will  be 
an  instantaneous  error  with  a  nonzero  variance. 

The  ocean  is  neither  homogeneous  nor  isotropic,  yet  measured  velocity 
turbulence  spectra  at  the  large  wave  numbers  tends  to  follow  the  power  law  relation 
that  is  applicable  to  homogeneous  and  isotropic  turbulence. 

E(k)  -  1.44  c2/3  k-5/3  (18) 

where  e  is  the  kinetic  energy  dissipation  rate  and  k  is  the  wave  number  of  the 
turbulent  eddies.  From  measured  data,  the  k"5/3  seems  to  hold  for  a  given  depth 
for  frequencies  as  small  as  1  cycle  per  10  hours. 

Assuming,  for  a  given  depth,  that  the  horizontal  velocity  components  satisfy 
the  conditions  for  equation  (18),  then  a  velocity  structure  function,  D(r),  can  be 
defined 

D(r)  =  fV (rj)  -  V(r2)]2  (19) 

where  V(r^)  and  V(r2)  are  the  instantaneous  velocity  components  at  positions  r^  and 
r2  and  where  r  =  |r^  -  r2|.  The  velocity  structure  function  is  related  to 
the  separation  distance,  r,  by  the  equation 

D(r)  =  1.  89  e2/3  r2/3  .  (20) 

Taking  the  square  root  of  equation  (20)  and  changing  units  of  r  to  nautical  miles  and 
\D(r)  to  knots  give 


Jb( r)  =  AV  =  1.52  e1/3  r1/3 


(21) 


♦The  value  of  0. 5  kn  for  the  maximum  current  due  to  internal  waves  is  selected  for 
illustration  purposes  and  is  probably  toward  the  upper  end  of  the  range  of  values 
encountered. 
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where 


AV  =  rms  velocity  differnece,  kn 
e  =  Energy  dissipation  rate,  cm2/s3 
r  =  Separation  distance,  nmi. 

In  the  oceans,  e  ranges  from  10“®  cm2/s3  to  10"!  cm2/s3  with  10~3  cm2/s3  being 
a  "typical"  or  representative  value.  Figure  6  gives  curves  for  these  three  values 
of  e  showing  how  the  rms  difference  between  the  readings  of  two  ideal  speed  sensors* 
separated  by  a  distance  r  would  vary. 

UNCERTAINTIES  CAUSED  BY  VERTICAL  CURRENT  GRADIENTS 

The  magnitudes  of  ocean  currents  characteristically  decrease  with  depth; 
surface  ships  and  submarines  at  various  depths  should  not  expect  to  experience  the 
same  local  water  movement.  (Current  directions  also  change  with  depth. ) 

Selected  current  velocity  profiles  appear  in  figure  7  for  three  stations  in  the 
Atlantic  Ocean  (selected  from  79  current  profiles  given  in  the  atlas2).  The  solid 
lines  represent  average  values  for  the  station  and  the  dashed  lines  represent  maximum 
values.  In  the  presence  of  gradients  such  as  those  shown  in  figure  7,  signif  icant 
differences  can  arise  in  measured  speed  relative  to  water  even  between  ships  of 
different  drafts,  say  a  carrier  and  its  shallow  draft  escorts.  A  submarine  operating 
at  a  depth  different  from  its  target  or  a  submarine  changing  depth  will  experience 
measured  speed  differences  caused  £  olely  by  the  current  gradients. 

The  histogram  in  figure  8  shows  the  difference  between  the  average  values 
of  surface  current  and  current  at  50  meters  depth.  Based  on  the  58  average  profiles 
in  the  North  Atlantic,  one  can  expect  about  0. 4  kn  average  difference  between  the 
magnitudes  of  surface  currents  and  currents  at  50  meters  depth.  Of  course,  for 
certain  specific  situations,  the  difference  could  be  as  much  as  3  to  5  knots. 

DISCUSSION 

The  local  water  movements,  characteristic  of  the  ocean,  will  cause  even  a 
well  designed  and  calibrated  water  speed  sensor  to  indicate  speed  readings  that  do 
not  accurately  represent  ship  motion.  For  example,  a  large  ship  maintaining  a 


*An  ideal  water  speed  sensor  is  one  that  senses  the  true  local  velocity  at  a  point 
without  perturbing  the  local  flow.  A  ship's  sensor  cannot  do  this  because  of  its 
finite  size  and  the  presence  of  the  ship's  hull. 
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constant  revolutions  per  minute  will  hold  a  relatively  steady  speed  over  ground.  *  The 
water  speed  sensor  measuring  the  local  water  velocity  at  the  sensor  will  measure 
this  steady  speed  plus  the  fluctuating  velocities  because  of  the  orbital  velocities  in 
the  surface  waves.  Given  sufficient  time,  these  indicated  speed  perturbations  ;an 
be  reduced  to  near  zero  by  appropriate  data  processing  and  averaging  techniques. 
However,  if  the  ship  is  maneuvering,  the  indicated  speed  must  follow  the  changes 
in  ship's  speed  to  be  useful.  The  time  and  space  scales  of  these  natural  oceanic 
water  movements  overlap  the  characteristic  time  and  space  scales  involved  in  many 
ship  maneuvers,  such  as  a  turn  taking  several  minutes  and  covering  a  mile  or  so  of 
path.  The  difficulty,  then,  is  how  to  obtain  speed  indications  that  represent  ship 
movements  rather  than  the  local  water  movements. 

The  magnitude,  frequency,  and  phase  of  the  difference  between  speed  relative 
to  water  and  speed  over  ground  are  dependent  on  the  vehicle  parameters,  depth  of 
the  sensing  region  below  the  surface,  geographical  location,  and  previous  weather 
conditions.  In  addition,  the  difference  between  the  "average"  water  speed  and  the 
instantaneous  water  speed  depends  on  the  "length"  of  the  averaging  period  or  how  far 
one  moves  through  the  "local  water  mass.  "  As  an  example  of  spatial  effects,  consider 
an  acoustic  doppler  speed  sensor  that  may  scatter  energy  from  water  several 
hundred  feet  deep  and  thus  measure  the  ship's  speed  relative  to  this  water.  Obviously 
the  effects  of  surface  waves  are  reduced,  but  then  one  must  contend  with  the  vertical 
current  gradients  (i.e. ,  the  change  of  magnitude  and  direction  of  local  ocean  currents 
with  depth).  The  measurement  and  use  of  "speed  relative  to  the  local  water  mass" 
requires  very  careful  consideration  of  what  the  "local  water  mass"  actually  represents. 

CONCLUSIONS 

•  Surface  waves  cause  oscillating  horizontal  currents  with  peak  magni¬ 
tudes  as  high  as  several  knots.  The  magnitudes  decrease  with  depth  with  the  spectral 
peaks  tending  toward  lower  frequencies  (or  longer  wavelengths).  This  spectral  change 
may  cause  difficulties  for  submarines  in  selecting  a  data  sampling  and  processing  or 
averaging  scheme  that  can  reduce  wave  effects  and  still  respond  to  ship  maneuvers. 

•  The  horizontal  currents  associated  with  internal  waves  are  propor¬ 
tional  to  the  wave  amplitude  and  phase  speed  and  have  magnitudes  up  to  tenths  of  knots. 
Since  the  wave  amplitudes  do  not  exceed  the  upper  layer  thickness  (depth  to  themocline) 
and  the  wave  phase  speed  can  be  calculated  or  measured,  it  should  be  possible  to 
estimate  the  upper  limit  to  the  horizontal  currents  associated  with  internal  waves  by 
using  local  measurements. 


*If  the  ship  changes  its  heading  relative  to  the  predominant  wave  field,  the  speed 
over  ground  may  change  by  several  percent  for  the  same  revolution  per  minute 
because  of  the  changing  drag  characteristics. 
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•  For  vessels  of  large  mass  In  relatively  rapid  local  water  movements, 
the  water  motion  causes  fluctuations  in  the  indicated  water  speed  while  the  speed  over 
ground  is  essentially  constant.  For  vessels  of  lesser  mass  in  slowly  varying  local 
water  movements,  the  speed  relative  to  water  is  very  nearly  constant  while  the  speed 
over  ground  varies  as  the  water  motion.  These  results  follow  directly  from  the 
inertial  and  drag  characteristics  of  the  vessel. 

•  Turbulence  causes  random  but  small  (~0. 01  to  0. 1  kn)  fluctuations 
in  the  speed  readout.  The  spectral  form  appropriate  to  homogeneous,  isotropic 
turbulence  may  be  useful  in  predicting  speed  readout  variations. 

•  Local  ocean  currents  tend  to  decrease  in  magnitude  and  shift 
directions  as  a  function  of  depth.  A  submersible  and  a  surface  support  vessel 
operating  together  will  almost  certainly  experience  differences  between  the  indicated 
water  speeds  of  the  two  vessels  even  while  they  maintain  zero  relative  velocity. 
Typical  gradients  for  the  North  Atlantic  show  about  0.4-kn  change  from  the  surface 
to  a  depth  of  50  meters. 
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Figure  1 

Velocity  Variations  Recorded  on  Board  Operational  Submarine  in  Sea  State  6  at 
Given  Depths  and  Speeds  (Signal  Bandpass  Approximately  0.  08  to  5  Hz) 
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- CALCULATED  CURVES  BASED  ON  EQUATION  (6) 

®  EXPERIMENTAL  VALUES  MEASURED  IN  ESTIMATED  SEA  STATE  6 


Figure  2 

RMS  Horizontal  Velocity  Caused  by  Wind  Generated  Surface  Waves 
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Figure  3 

Examples  of  Density  Gradients  in  Atlantic  Ocean3 
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Figure  5 

Variations  of  Speed  of  Submersible  in  Oscillatory  Current  Due  to  Internal  Waves 
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Figure  5  (Cont) 
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RMS  VELOCITY  DIFFERENCE,  kn 


Figure  6 

RMS  Difference  Between  Water  Velocity  at  Two  Points  Separated  by  a  Distance  r 
in  an  Isotropic,  Homogeneous  Turbulence  Field  Characterized  by  Dissipation  Rate,  e. 
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Figure  8 

Histogram  of  Difference  between  Current  Magnitude  at  Surface  and  50-meter 
Depth  for  58  Stations  in  Atlantic  Ocean. 
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APPENDIX  A 


SPEED-MEASUREMENT  EFFECTS  CAUSED  BY 
OSCILLATORY  CURRENTS 

The  oscillatory  currents  in  an  internal  wave-field  will  cause  a  variation  of 
indication  from  a  speed  sensor  that  is  measuring  vehicle  speed  relative  to  the  water. 
The  drag  forces  on  the  submersible  vehicle  depend  on  its  velocity  relative  to  the 
water,  but  the  acceleration  forces  depend  on  velocity  changes  relative  to  the  ground. 
The  oscillating  horizontal  currents  of  internal  waves  affect  both  speed  relative  to 
water  and  speed  over  the  ground. 

Consider  a  submersible  of  mass  M  maintaining  a  constant  thrust  force  TQ  and 
moving  in  a  region  of  oscillatory  horizontal  currents  v  given  by 

v  =  v0  sin  wt 


where: 


v0  =  the  peak  current  value 

w  =  the  encounter  frequency  seen  by  an  observer  on  the  submersible 
t  =  time. 

Assuming  the  drag  force  is  proportional  to  the  square  of  the  speed  relative  to  the 
water,  the  equation  of  motion  is 

M  ~  +  Cd(V-v)2  +  T0  =  0 


where: 

M  =  mass  of  submersible 
V  =  speed  of  submersible  relative  to  ground 
Cd  =  drag  force  coefficient 
v  -  current  relative  to  ground 
Tq  -  constant  thrust  force. 


_  Zo  - 
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To  obtain  a  solution,  note  that  if  v  =  0,  the  constant  thrust  TQ  balances  the  drag 
at  steady  speed  V0  such  that 

Cd  V02  +  T0  =  0  . 

The  speed  over  ground  V  can  be  expressed  as  the  steady  speed  VQ  plus  a  perturbation  u, 

V  =  V0  +  u  . 

Substituting  this  into  the  governing  differential  equation  and  rearranging  ter  ns  gives 

M^+Cd  [2V0(u-v)  +  (u-v)2]  =  0  . 

Assuming  u-v  «  V0,  the  equation  simplifies  to 

+2CdV0(u-v)  =  0  . 

Substituting  for  v  and  rearranging  now  gives 


dt  M  M  ° 


The  solution  to  this  equation  is 


u  =  vQ  cos  <p  sin  <wt  +  tp) 


where 


tan<P  =  - 


2V0Cd  * 


The  speed  over  the  ground  is  given  by  V  =  V0  +  u  or 
V  =  VQ  +  vQ  cos  V  sin  (wt  +  <P)  . 

The  indicated  speed  (or  speed  relative  to  the  water)  is  given  by  Vj  =  V  -  v  or 


Vj  =  VQ  +  v0  sin  <Pcos(u>t  +  <P)  . 
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